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herein  has  been  reported  in  BRL  Reports  and  Memorandum  Reports.  The 
thrust  of  this  present  effort  has  been  the  extraction  and  discussion 
of  ionospheric  effects  on  radio  propagation  which  continue  to  be  of 
importance  in  the  performance  and  effectiveness  of  communications  sys- 
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I INTRODUCTION 


In  the  years  1948-71,  the  Ballistic  Research  Laboratories  (BRL) 
developed  several  methods  for  the  study  of  ionospheric  characteristics 
with  the  principal  objective  of  understanding  the  effects  of  the  ionosphere, 
disturbed  and  undisturbed,  upon  the  propagation  of  radio  signals  used  for 
the  tracking  of  objects  emitting  signals  and  in  the  transmission  of  in- 
formation. The  major  efforts  focused  on  the  use  of  radio  methods  for  the 
diagnostic  studies  in  the  belief  that  the  most  direct  access  to  propagation 
effects  could  thereby  be  obtained.  This  approach  of  necessity  led  to  in- 
ferential determination  of  ionospheric  characteristics  through  the  postu- 
lation of  simple  or  complex  models  of  ionospheric  structure  followed  by 
iterative  modifications  of  the  models  to  provide  compatibility  between 
predicted  and  observed  propagation  effects.  Sensitivity  of  the  method  was 
provided  through  use  of  multiple  coherent  frequencies,  careful  attention 
to  antenna  and  polarization  properties  of  the  equipments  and  sophisticated 
detection  and  recording  of  the  received  signals. 

The  earliest  efforts  of  the  BRL  to  assess  the  effects  of  the  ionosphere 
on  radio  signals  were  associated  with  the  mission  of  developing,  installing 
and  operating  the  missile  tracking  system  known  as  Doppler  Velocity  and 
Position  (DOVAP)  at  the  White  Sands  Missile  Range  in  the  late  1940s.  This 
system,  which  will  be  described  briefly  in  a following  section,  was  capable 
of  real-world  tracking  accuracies  equivalent  to  the  best  of  the  ballistic 
cameras  and  was  not  equaled  in  accuracy  by  microwave  systems  until  the  ad- 
vent of  sophisticated  phased-array  radars  in  the  late  1960s.  A precursor 
to  the  DOVAP  system  was  developed  in  Germany  in  the  early  1940s  to 
provide  accurate  velocity  information  for  long  range  ballistic  rockets  (V-2) 
since  radar  in  those  years  was  incapable  of  providing  the  needed  data. 

With  the  DOVAP  system,  basic  trajectory  information  is  derived  by  measur- 
ing, at  the  ground,  the  phase  of  radio  signals  transmitted  to 
and  from  the  object  being  tracked.  Distances,  then,  are  related  to 
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phase  through  the  transmitted  frequency  and  velocity  of  propagation  (phase 
velocity) . The  phase  velocity  of  propagation  in  the  ionosphere  is  greater 
than  the  vacuum  velocity  of  light  so  that  DOVAP  range  information  is  fore- 
shortened with  respect  to  true  values  for  objects  tracked  in  and  above  the 
ionosphere.  For  the  DOVAP  system  operating  at  White  Sands  Missile  Range 
with  frequencies  of  37  and  74  MHz,  the  range  error  could  be  as  much  as  2 
kilometers  at  an  altitude  of  1,000  kilometers.  Efforts  at  the  BRL  were 
directed  to  finding  ionosphere  corrections  for  the  DOVAP  data;  by  a modeling 
process,  these  corrections  could  be  used  to  describe  the  altitude  profile  of 
total  electron  content  in  the  ionosphere.1  The  same  basic  method,  but  with 
much  improved  data  handling  and  analysis,  was  developed  for  determining  the 
ionospheric  structures  of  Mars  and  Venus  during  the  Mariner  missions.2  J 4 5 

In  the  late  1950s,  the  two-frequency  method  of  determining  ionospheric 
properties,  developed  first  by  Seddon,e  was  adopted  by  the  BRL  to  enable 
more  accurate  measurements  of  ionospheric  properties  and  to  free  such 
measurements  from  restriction  to  launch  facilities  equipped  with  the  DOVAP 
tracking  system.  The  BRL  utilized  frequencies  of  36/145  MHz  for  the  two- 
frequency  experiment  because  of  long  familiarity  with  appropriate  equipments, 
because  data  from  frequencies  below  approximately  20  MHz  become  increasingly 
difficult  to  interpret,  because  rocket  antennas  are  relatively  easy  to  design 
for  such  frequencies,  and  because  of  minimum  background  interference. 

This  adopted  method  was  applied  in  a series  of  experiments,  identified 
as  STRONGARM,  carried  out  at  Wallops  Island,  Va. , in  1959-60. 7 Shortly  after 
these  experiments,  a crash  effort  was  initiated  to  prepare  and  field  rocket- 
borne  multi-frequency  beacons  for  the  high  altitude  nuclear  tests  at  Johnston 
Island  in  1962. 8 For  this  program,  a third  coherent  frequency  at  888  MHz 
was  added  to  provide  greater  dynamic  range  in  the  measurements  taken. 

Eighteen  rockets  equipped  with  these  multi-frequency  beacons  were  flown 
during  the  test  program. 


* The  frequencies  noted  throughout  this  report  are  rounded  to  the  nearest 
MHz.  Where  frequency  ratios  are  given,  such  frequencies  are  phase-coherent. 


Following  the  nuclear  tests  at  Johnston  Island,  the  BRL  participated  in 
the  Test  Readiness  Program  called  for  under  Safeguard  C of  the  Limited 
Test  Ban  Treaty.  Supported  by  this  program,  the  multi-frequency  propagation 
experiment  was  greatly  improved  in  ground  instrumentation,  data  processing 
and  analysis,  antenna  design,  and  rocket-borne  transmitting  equipment. 

Eventually,  five  coherent  frequencies  were  incorporated  into  the  experiment 
design  covering  the  range  36-875  MHz.  Under  the  Test  Readiness  Program,  the 
BRL  propagation  experiment  was  brought  to  a high  degree  of  perfection. 

In  the  years  after  1962,  ionospheric  measurements  by  the  BRL  were 
always  carried  out  in  support  of  field  programs  whose  objectives  were  sub- 
stantially broader  than  the  isolated  measurement  of  ionospheric  parameters. 

Thus,  the  BRL  provided  ionospheric  data  for  the  barium  release  programs 
SECEDE  II  and  SECEDE  III9  10  and  the  coordinated  geophysical  studies  in 
the  polar  cap  absorption  (PCA)  event  in  November  1969. 11  12  All  the  high 
altitude  rocket  experiments  carried  out  by  the  BRL  have  been  catalogued;13 
however,  the  more  significant  data  samples  of  ionospheric  electron  content 
and  densities  from  rocket  experiments  are  listed  in  Tables  1.1,  1.2,  and  1.3. 

In  time,  data  taken  in  the  auroral  regions  at  times  of  sunspot  maximum 
(including  the  data  in  Tables  1.1  and  1.2)  may  become  increasingly  important 
if  sunspot  activity  is  entering  a declining  phase.  Certainly  solar  activity 
in  the  1968-69  period  }f  sunspot  maximum  was  considerably  less  than  the 
1957-58  period,  and  it  has  been  noted  recently  that  the  11-year  solar 
cycle  may  not  be  a constant  characteristic  of  the  sun  and  solar  activity  has  been 
very  small  for  long  periods  of  time  in  the  past.14  Thus,  design  of  communica- 
tions systems  for  the  relatively  benign  radio  transmission  environment 
accompanying  times  of  low  sunspot  activity  could  yield  marginal  or  disrupter 
performance  in  the  more  vigorous  ionospheric  environment  created  by  solar 
or  man-made  disturbances  in  the  high  atmosphere. 

In  the  late  1950s,  the  BRL,  as  well  as  others,  became  interested  in 
the  use  of  satellites  as  carriers  for  radio  transmitters  and  the  promise 
of  a major  new  capability  for  ionospheric  studies  using  signals  from  such 
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Data  Type:  [1]  Total  electron  Content,  Trajectory  Comparison;  [2]  Total  Content,  Faraday  Rotation;  [J]  Interred 

Electron  Densities. 

Altitude  interval  (kilometers)  witliin  wliieh  data  were  obtained. 
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transmitters.  When  the  Sputnik  I was  launched,  the  BRL  had  designed  and 
was  already  assembling  ground  receiving  equipments  for  use  with  satellite 
transmitters  and  was  quickly  able  to  initiate  tracking  of  Sputnik.  This 
early  experience  made  possible  a rapid  development  of  data  reduction  methods 
and  analysis,  and  incorporation  of  the  necessary  receivers  and  antennas 
in  the  ground  station  used  for  the  multi-frequency  rocket  experiments.  Thus, 
ionospheric  measurements  could  be  made  as  soon  as  the  two-frequency  navy 
navigation  satellites,  launched  under  the  Transit  Program,  were  placed  in 
orbit.  When  the  U.S.  high  altitude  nuclear  tests  were  carried  out  in  the 
Pacific  in  1962,  it  was  possible  to  measure  total  electron  content  in  the 
ionosphere  (below  approximately  900  kilometers)  for  95  passes  of  Transit 
4/\  and  Transit  4B  (Omicron  1 and  Alpha  Eta  1)  in  the  period  17  May  through 
4 November  1962.  The  ground  tracking  station  was  located  at  Johnston  Island 

V 

and  data  were  obtained  at  latitudes  of  5-25  degrees  norLh  and  longitudes 
of  160-180  degrees  west."  lj  The  methods  developed  at  the  BRL  and  by  other 
investigators  made  possible  ionosphere  corrections  to  position  and  velocity 
information  for  the  many  users  of  the  Transit  navigation  system. 


For  a variety  oi  reasons,  a major  one  being  changing  roles  and  missions, 
BRL  activities  in  rocket  measurements  in  the  upper  atmosphere  were  phased 
out  rather  quickly  after  1971.  The  major  contribution  of  the  high  altitude 
rocket  measurements  and  satellite  tracking  programs  of  the  BRL  was  the 
development  and  perfection  of  methods  for  analyzing  and  understanding  the 
transmission  of  phase-coherent,  multi-frequency  signals  through  disturbed 
ionospheres . 


Pi 


II  METHODS  EMPLOYED 

As  mentioned  in  the  previous  section,  the  thrust  of  BRL  efforts  to 
determine  ionospheric  effects  on  radio  signals  changed  from  an  emphasis  on 
correcting  radio  tracking  data  to  the  use  of  radio  signals  for  diagnostic 
measurements  of  disturbed  ionospheres.  This  changing  emphasis,  together 
with  rapidly  improving  electronic  and  recording  equipments,  led  to  evolu- 
tionary changes  in  instrumentation  as  well.  The  result  was  a better  capa- 
bility to  assess  ionospheric  effects  and  provide  a base  of  information  to 
predict  these  effects  for  a variety  of  ionospheric  structures  and  other 
frequency  regimes. 

In  this  section,  the  methods  used  by  the  BRL  in  assessing  ionospheric 
effects  on  radio  signals  are  briefly  reviewed;  more  detailed  discussions  are 
found  in  various  BRL  reports.15  16  17  1S 

A.  Trajectory  Comparison  Method 

As  mentioned  in  the  previous  section,  the  earliest  BRL  efforts  in 
investigation  of  ionospheric  effects  were  associated  with  the  operation  of 

| 

the  DOVAP  missile  tracking  system  at  the  White  Sands  Missile  range.  The 
^ ionospheric  effects  associated  with  the  phase  comparison  tracking  system  were 

later  used  for  inferential  determination  of  ionospheric  structures  of  Mars 
and  Venus.2  3 5 

A portion  of  the  DOVAP  tracking  system  is  represented  schematically 
in  Figure  2.1.  In  operation,  a radio  signal  from  the  transmitter  is  received 
by  the  missile,  doubled  in  frequency  and  retransmitted  to  three  or  more 
ground  receivers.  The  same  transmitter  signal  is  sent  by  ground  link  to  the 
same  receiver  sites,  doubled  in  frequency  and  mixed  with  the  signals  re- 
ceived via  the  rocket  transmitter.  Because  phase  coherence  is  retained 
throughout  the  system,  the  resultant  total  signal  phase,  after  mixing  and 
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Figure  2.1  SCHEMATIC  OF  RADIO  TRACKING  SYSTEM 


multiplication  by  the  wavelength  of  the  doubled  frequency,  provides  a 
measure  of  the  distance  from  the  transmitter-to  rocket-to  ground  receiver. 
Knowledge  of  the  total  phase  requires  continuous  recording  of  the  mixed 
signals  from  the  time  of  missile  lift-off.  For  the  general  case  of  non- 
coincident transmitter  and  receiver,  the  measured  distance  generates  an 
ellipsoid  of  revolution.  The  intersection  of  three  such  ellipsoids  (re- 
quiring three  receivers)  locates  the  missile  in  space.  For  the  special  case 
of  coincident  transmitter  and  receiver,  the  ellipsoid  becomes  a sphere. 

Because  rocket-borne  antennas  frequently  have  dipole  patterns,  ground 
antennas  have  circular,  or  quasi-circular,  polarization  to  avoid  signal 
dropouts  and  abrupt  phase  shifts.  As  a consequence,  the  phase  of  received 
signals  is  affected  by  missile  spin  and  polarization  rotation  arising  from 
the  Faraday  effect  (earth's  magnetic  field  interaction  with  the  electric  field 
of  the  radio  signals).  Accordingly,  reduction  of  measured  phase  accumulation 
to  determine  missile  position  must  consider  all  phase  changes  arising  from 
t-vector  rotation. 


In  the  reduction  of  phase  data  acquired  by  the  DOVAP  system,  the 
assumption  is  made  that  the  wavelength  used  in  converting  phase  to  distance 
is  that  appropriate  to  signal  propagation  in  vacuo,  ie.,  Af=Co  when  and  f 
are  the  wavelength  and  frequency  of  the  tracking  signal  used,  and  is  the 
vacuum  speed  of  light.  In  fact,  however,  in  a dispersive  medium  such  as 
the  ionosphere,  the  propagation  speed  is  not  C^;  for  a phase-measuring 
system  the  phase  velocity  of  propagation  in  the  ionosphere  is  greater  than 
Cq.  Thus,  the  range  measurements,  S^,  derived  in  the  DOVAP  system  are  given 
by  the  relation,  S^=4>pCQ/f  (where  is  the  measured  total  phase  at  any  point 
in  the  trajectory)  when  in  fact  the  true  range,  S,  is  given  by,  S=;^c/f, 
with  C being  the  correct  value  for  propagation  speed  in  the  dispersive 
medium.  The  relationship  of  C to  the  various  parameters  describing  the  iono- 
sphere (electron  and  ion  densities,  collision  frequencies,  electron  and  ion 
gyrof requencies , etc.)  has  been  known  for  many  years  and  discussed  ex- 
tensively.19 If  the  basic  system  frequency  is  sufficiently  high,  the  true 
speed  of  propagation  is  almost  entirely  a function  of  the  free  electron 


density  in  the  ionosphere  and  can  be  expressed  by  a very  simple  relationship. 
Farther,  for  sufficiently  high  frequencies,  the  assumption  of  straight- 
line  propagation  in  the  ionosphere  introduces  negligible  error  in  derived 
trajectory  parameters.  If  the  ionosphere  has  horizontal  uniformity  in  the 
space  occupied  by  the  missile  and  the  tracking  system,  differences  in 
and  S can  be  related  simply  to  altitude  differences,  h^  and  h.  Thus,  a 
measurement  of  the  quantities  and  S,  together  with  a simple  expression 
relating  their  difference  to  ionospheric  electron  density,  makes  possible, 
through  an  inversion  process,  the  measurement  of  electron  content  in  the 
ionosphere  below  the  altitude  of  a tracked  missile. 

In  the  absence  of  an  independent  range  measuring  system,  unaffected  by 
the  ionosphere,  the  quantities  S or  h are  not  directly  available.  However, 
the  DOVAP  system  is  sufficiently  accurate  that  a vacuum  trajectory  (or 
orbit)  may  be  computed  using  initial  conditions  from  t lie  DOVAP  trajectory 
after  missile  burnout  and  above  the  effective  drag  region  but  before  the 
missile  enters  the  ionosphere.  The  quantities  S or  h are  given  by  tiie  vacuum 
trajectory  parameters.  In  Figure  2.2  is  shown  a typical  profile  of  electron 
density  versus  altitude  derived  through  application  of  the  method  described. 

The  error  bars  are  representative  of  the  accuracy  of  this  method.  Accuracies 
are  poor  for  small  electron  densities  and  small  missile  speeds. 

B.  Faraday  Rotation 

Because  of  the  earth's  magnetic  field,  the  ionosphere  is  birefringent 
so  that  a linearly  polarized  wave  propagating  through  the  ionosphere  is 
broken  into  two  counter-rotating  circularly  polarized  waves  traveling  at  slightly 
different  velocities.  These  are  familiarly  known  as  the  ordinary  and  extra- 
ordinary waves.  The  combination  of  these  two  waves,  as  seen  by  an  antenna, 
as  the  waves  propagate  would  appear  to  be  a linearly  polarized  wave  whose  plane 
of  polarization  is  slowly  rotating.  This  rotation  is  called  Faraday  rotation 
after  its  discoverer.  The  same  phenomenon  occurs  for  circularly  polarized 
waves  as  well.  In  effect,  the  Faradiy  rotation  advances  or  retards  the 
rotating  vector  as  the  wave  proceeds.1' 
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As  mentioned  earlier,  missile  spin  and  Faraday  rotation  introduce 
phase  contributions  to  the  total  measured  phase,  c , which  must  be  removed 
in  providing  accurate  trajectory  information.  To  provide  the  needed  cor- 
rections, each  receiver  site  acquires  phase  data  through  two  oppositely 
circularly  polarized  antennas.  Mixing  of  the  signals  received  through 
the  two  antennas  provides  phase  data  proportional  to  polarization  rotation 
resulting  from  the  sum  of  missile  spin  and  Faraday  rotation.  This  phase 
information  permits  the  necessary  corrections  for  accurate  trajectory  deter- 
mination. Through  a simple  analytical  procedure  it  is  possible  to  separate 
the  missile  spin  and  Faraday  contributions  to  the  total  E-vector  rotation.18 
In  fashion  analogous  to  the  trajectory  comparison  technique,  missile  spin 
is  determined  above  the  altitude  of  appreciable  drag  but  prior  to  missile 
entry  into  the  ionosphere.  Figure  2.3  illustrates  the  method  of  deter- 
mining Faraday  rotation  information  from  the  recorded  DOVAP  spin  data. 

For  frequencies  used  in  the  trajectory  comparison  method,  Faraday  rotation 
is  proportional  to  the  total  electron  content  between  the  missile  and 
ground  receivers.  Determination  of  electron  content  from  Faraday  rotation 
data  is  less  accurate  than  information  from  trajectory  comparison;  however, 
the  data  are  relatively  easy  to  obtain. 

C . Dispersive  Doppler 

Referring  to  Figure  2.1,  it  is  clear  that  if  the  information  desired 
is  the  structure  of  the  ionosphere,  the  transmission  of  two  or  more  fre- 
quency-coherent signals  from  the  missile  will  provide  data  more  simply  and 
accurately  than  the  trajectory  comparison  method.  Choosing  a relatively 
large  frequency  ratio  (although  this  is  certainly  not  necessary)  and  the 
lowest  frequency  such  that  differences  in  ray  paths  can  be  ignored,  disper- 
sive or  differential  phase  differences  are  readily  obtained  by  mixing  the 
signals  received  at  the  ground,  after  suitable  multiplication  of  the 
lower  frequency  by  the  frequency  ratio.  This  dispersive  phase  can  be 
readily  related  to  ionospheric  parameters  and  the  techniques  have  been 
extensively  discussed.6  16  17  20  All  BRL  activities  after  1958  related  to 
ionospheric  propagation  measurements  employed  the  dispersive  doppler  or 
differential  phase  technique. 
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Ionospheric  measurements  using  the  techniques  of  Faraday  rotation  and 
dispersive  doppler  have  been  obtained  with  signals  transmitted  by  satellites. 

For  satellites,  it  is  not  easily  possible  to  maintain  continuous  signal 
reception  from  lift-off  of  the  satellite  so  that  a phase  ambiguity  of  many 
cycles  exists  in  the  dispersive  doppler  or  Faraday  rotation  data  for  any 
given  pass  of  the  satellite  over  a ground  receiving  station.  For  the  case 
of  Faraday  rotation,  the  ambiguity  may  be  removed  m theor>  for  a given 
satellite  pass  if  signal  propagation  normal  to  the  earth's  magnetic  field 
occurs  sometime  during  that  pass.  For  such  a condition,  the  Faraday  rotation 
is  zero  so  that  accumulated  Faraday  cycles  can  be  determined  by  counting 
forward  and  backward  from  the  normality  condition.  In  practice,  however, 
the  propagation  vector  is  not  normal  to  the  earth's  field  along  the  entire 
propagation  path  so  at  best  a rather  fuzzy  region  occurs  in  a portion  ol 
track  where  the  ray  path  is  close  to  the  perpendicular  configuration;  and 
the  ambiguity  is  not  easily  resolved.  If,  however,  dispersive  doppler  and 
Faraday  rotation  are  recorded  simultaneously,  it  is  possible,  by  a hybrid 
method,  to  remove  ttie  ambiguity  and  derive  data  for  total  electron  content 
below  the  satellite  altitude.15  1 The  hybrid  technique  was  used  to 

develop  an  extensive  data  set  during  the  high  altitude  nuclear  tests  in 
1962,  referred  to  earlier  in  Section  I (and  see  Figure  3.1). 

) . Signal  At tenuation 

As  a general  statement,  signal  attenuation  derives  from  several  dis- 
tinct mechanisms.  The  familiar  and  rather  obvious  mechanism  is  that  of 
spherical  expansion  of  a transmitted  wave,  with  the  energy,  or  power,  de- 
creasing as  tiu  inverse  square  of  the  distance  from  the  transmitter. 

linos t as  familiar  as  spherical  expansion  is  the  absorption  of  signal 
energy  through  tin  collision  of  electrons  oi  ions  witli  each  other  or  with 
jtral  const! tu<  t in  the  atmosphere  luring  tassagc  ol  th<  signal. 

Tin  equations  whim  govern  absorption  ire  well  known.19  Thi  ••  docs  not  imply 
that  tin.  atmospheric  parameters  are  sufficiently  w ii  known  that  absorp- 
tion can  be  confidently  predicted  in  all  frequency  regimes  or  for  all  levels  * 

o£  atmospheric  disturbance.  The  third  source  oi  attenuation  arises  from 

r 
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a variety  of  propagation  conditions  which  may  amplitude-  and  phase-modulate 
the  passing  signals  (e.g.,  diffraction  and  multipath)  or  for  which  the  radius 
of  curvature  of  an  expanding  wave  front  is  markedly  changed  (e.g.,  focusing 
and  defocusing) . This  last  source  of  attenuation  will  be  discussed  more 
fully  in  the  following  sections. 

Signal  attenuation  arising  from  the  absorption  of  energy  is  very  impor- 
tant for  significant  portions  of  the  radio  spectrum  (HF  frequencies  and 
below  in  the  ionized  atmosphere,  X-band  and  above  in  the  non-ionized  at- 
mosphere). The  method  employed  by  the  BRL  was  quite  s traight-forward. 

The  zero-absorption  level  was  defined  as  the  signal  level  (properly  cor- 
rected for  antenna  pattern  effects)  at  the  points  of  entry  into  and  exit 
from  the  ionosphere.  A straight  line  joining  these  points,  in  a plot  of 
signal  level  (db  below  one  milliwatt)  as  a function  of  slant  range,  es- 
tablished the  zero  absorption  signal  level  during  the  time  the  tracked 
vehicle  was  in  or  above  the  ionosphere.  Signal  level  departures  from  this 
straight  line  interpolation  were  attributed  to  absorption  of  the  signal 
energy.  An  example  of  measured  absorption  derived  by  this  method  is  shown 
in  Figure  3.2,  where  slant  range  has  been  converted  to  altitude  of  the 
vehicle  being  tracked. 


20 


Ill  PROPAGATION  EFFECTS 

In  the  previous  sections,  reference  was  made  to  effects  of  the  ionosphere 
on  a particular  radio  tracking  system  (DOVAP)  used  for  trajectory  determina- 
tions. In  point  of  fact  there  are  ionospheric  effects  on  such  systems 
which  were  not  described;  such  effects  are  important  for  a much  broader 
family  of  radio  systems  used  for  communications  purposes.  Further,  iono- 
spheric effects  are  important  for  radio  navigation  systems  and  the  greater 
the  demand  for  system  accuracy,  the  more  important  luvones  the  knowledge 
of  propagation  factors  required  for  ioiv  eiuri  i n.  flu  basic 

methods  described  in  Section  II  are  thos.  •.  < t ....  at  e it!)  major  improve- 
ments in  precision,  accuracy  and  data  hand1.  . . . this  section, 

several  propagation  effects  not  addressed  ■ 1:  previous 

sections  will  be  discussed. 

A.  Ionospheric  Tilts 


For  purposes  of  determining  propagation  effects  in  the  various  radio- 
frequency systems,  the  ionosphere  is  frequently  assumed  to  be  horizontally 
stratified.  Although  such  an  assumption  simplifies  calculations  and  ex- 
pressions for  propagation  path  refraction,  marked  horizontal  gradients  in 
ionization  occur  on  occasion  everywhere  on  earth  and  very  frequently  for 
a large  portion  of  the  ionosphere.  The  so-called  ionospheric  "tilts"  make 
difficult  the  prediction  of  usable  frequencies  for  given  point-to-point 
communications  at  high  frequency  (HF) , are  responsible  for  anomalous  HF 
propagation,  and  give  rise  to  the  phase  and  amplitude  variations  of  impor- 
tance for  VLF  and  ELF  communications  systems.  In  recent  years,  the  burden 
of  radio  communications  has  shifted  to  higher  capacity,  higher  frequency 
systems,  but  abundant  interest  remains  in  high  frequency  communications 
on  the  part  of  many  present  and  potential  subscribers.  Ionospheric  tilts 
can  enhance  as  well  as  degrade  the  range  of  HF  systems,  and  can  give  rise 
to  the  irritating  condition  of  non-reciprocal  communications. 
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An  example  of  a marked  ionospheric  tilt  is  shown  in  Figure  3.1.  The 

data  were  obtained  in  1962  from  the  tracking  of  Transit  satellites  using 

the  hybrid  Faraday-Doppler  technique.""  A twofold  increase  in  total  electron 

content  occurred  between  the  north  latitudes  of  18°-22°,  or  west  longitudes 

of  163°-155°,  or  some  combination  of  both.  At  the  local  time  of  satellite 

passage,  gradients  in  the  west-east  direction  are  normally  small  and  the 

observed  gradient  is  probably  indicative  of  a strong  south-to-north  increase 

in  content.  In  the  summer  northern  hemisphere,  a maximum  in  electron  content 

usually  occurs  between  15°-25°,  arising,  most  likely , from  subsidence  of 

> > 

ionized  constituents  convected  to  high  altitudes  by  E x B drift  in  the 
equatorial  electrojet.  The  data  shown  in  Figure  3.1  indicate  a higher- 
than-usual  maximum  electron  content.  From  the  standpoint  of  HF  propagation, 
the  electron  content  profile  shown  would  indicate  a substantial  ionospheric 
tilt  and  a difficult  situation  for  establishing  a north-south  communications 
link. 


B.  Radio  Signal  Absorption 

In  the  previous  section,  sources  of  signal  attenuation  in  the  iono- 
sphere were  briefly  summarized.  The  absorption  of  signal  energy  occurs  under 
conditions  when  the  signal  frequencies  are  of  the  same  order  as  the  col- 
lision frequencies  of  the  ionized  particles  with  the  neutral  atmospheric 
particles,  and  with  each  other.  Absorption  also  occurs  when  the  signal 
frequency  approaches  the  plasma  frequency  of  the  ionized  medium.  For 
the  latter  case,  however,  ionization  gradients  become  very  important  because 
refraction  of  the  signals  becomes  very  large  and  the  signals  are  essentially- 
reflected  away  from  the  strongly  absorbing  regions.  Thus,  VLF  and  ELF  are 
trapped  in  the  so-called  ionospheric  waveguide  with  relatively  little  energy 
removed  by  electron  collisions. 
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For  HF,  and  for  lower  frequencies  as  well,  electron  collisions  can 
remove  substantial  energy,  and  absorption  is  an  important  consideration  w n 
appreciable  ionization  exists  at  altitudes  where  electron  collision  frequen- 
cies approach  signal  frequencies.  In  disturbed  ionospheres,  in  particulsr 
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the  aurorally  disturbed  ionosphere,  HF  signal  absorption  is  very  significant. 
In  a series  of  experiments  during  the  polar  cap  absorption  (PCA)  event  in 
November  1969,  the  ERL  measured  Hi’  and  VHF  absorption  at  intervals  during 
the  event.  One  such  record  is  shown  in  Figure  3.2.  Because  collision 
frequency  is  strongly  dependent  on  altitude,  observed  absorption  is  very 
sensitive  to  the  altitude  distribution  of  free  electrons  as  well  as  the 
frequency  employed.  Further,  because  of  the  earth's  magnetic  field,  the 
amount  of  absorption  depends  on  the  polarization  of  the  transmitted  signals. 

C.  Signal  Refraction 

As  mentioned  above,  signal  refraction,  or  bending,  becomes  marked  when 
signal  frequencies  approach  the  plasma  frequency  of  the  medium.  Under  cir- 
cumstances usually  encountered  in  propagation  problems,  refraction  occurs 
with  little  loss  in  signal  energy  (i.e.,  collision  frequencies  are  much  smalle 
than  the  signal  frequency) . Ionization  gradients  determine  the  radius  of 
curvature  of  the  refracted  signal.  Refraction  is  of  great  importance  at 
HF  and  lower  frequencies  since  it  provides  the  mechanism  for  long-distance 
point-to-point  communication  for  earth  terminals.  Refraction  is  of  importance 
also,  for  radio  systems  employing  frequencies  greater  than  HF  for  which 
line-of-sight  bearing  is  very  important.  Target  tracking  radars  fall 
into  this  important  class  of  systems.  Signal  refraction,  albeit  small,  pro- 
duces significant  target  position  errors.  A technique,  analogous  to  the  dis- 
persive doppler  technique  discussed  in  Section  II,  has  been  used  to  correct 
line-of-sight  errors  in  target  tracking  systems.  Two  frequencies,  with  suf- 
ficient frequency  separation,  are  used  in  the  line-of-sight  measurement. 
Measured  angle  differences  in  the  lines-of-sight , together  with  an  approxi- 
mate model  of  the  refracting  medium,  provide  a method  for  estimating  with 
good  accuracy  the  true  bearing,  or  line-of-sight,  of  a target  being  tracked. 

As  in  the  trajectory  comparison  method  described  in  Section  II,  the  use  of 
the  affected  system  to  determine  self-inflicted  ionospheric  effects  is  a 
useful  and  effective  method  for  a pragmatic  determination  and  subsequent 
correction  of  errors  caused  by  such  effects. 
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(SECONDS  AFTER  LAUNCH) 


D . Focusing,  Defocusing,  Multipath  and  Diffraction 

Propagation  problems  falling  under  this  general  classification  produce 
difficulties  in  all  frequency  regimes.  Ionospheric  discontinuities  or 
significant  gradients  of  size  comparable  to  a few  wavelengths  of  the  trans- 
mitted signals  produce  large  amplitude  and  phase  fluctuations  in  received 
signals.  Larger  discontinuities  can  produce  significant,  but  slower  varying 
amplitude  changes;  because  accompanying  phase  changes  are  also  slow, 
effects  on  phase  sensitive  systems  are  not  ..larked  unless  system  margins  are 
small  (the  order  of  a few  db). 

Ionospheric  perturbations  caused  by  rapidly  varying  particle  precipi- 
tation (auroral  regions  and  nuclear  explosions),  plasma-producing  chemical 
releases  such  as  barium,  and  nuclear-generated  plasmas  give  rise  to  the 
propagation  problems  discussed  in  this  section.  In  addition  to  measure- 
ments of  ionospheric  perturbations  caused  by  high-altitude  nuclear  explosions, 
the  BRL  has  carried  out  numerous  propagation  studies  in  the  auroral  regions  and 
associated  with  barium  releases.10  23  24  For  the  latter,  there  are  many  and 
clear  examples  of  the  propagation  effects  relevant  to  the  problems  here 
discussed.  Some,  but  not  all,  of  the  relevant  BRL  experiments  are  listed 
in  Tables  1.2  and  1.3. 

The  discussion  here  will  be  limited  to  illustration  of  the  propaga- 
tion problems  associated  with  ionized  barium  clouds.  In  Section  IV,  an 
illustration  of  auroral  disturbance  will  be  discussed.  In  Figure  3.3 
is  shown,  schematically,  the  effect  of  a plasma  "ball"  (approximating 
a barium  ion  cloud)  on  radio  signals  transmitted  in  close  proximity.  The 
region  indicated  as  "multipath"  is  also  a region  of  focusing  and  diffrac- 
tion. In  October  1967,  the  BRL  participated  in  an  Air  Force-sponsored  chemical 
release  program,  HAVEN  BARTER.  The  BRL  experiment  consisted  of  a three- 
frequency  phase-coherent  transmitter  carried  on  the  same  vehicle  from  which 
the  chemicals  (barium)  were  released.  The  transmitter  was  "viewed"  by  two 
ground  stations  such  that  signals  transmitted  through  and  around  the  several 
barium  ion  clouds  were  observed.  Figure  3.4  is  a record  of  the  72  MHz 
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Figure  3.4  SIGNAL  LEVEL  OF  72MHz  TRANSMISSION  RECEIVED  AT  WALLOPS  ISLAND  LAUNCH  SITE  PROJECT  HAVEN  BARTER 


; 
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signal  received  by  the  ground  station  covered  by  the  "footprint"  of  the 
ion  clouds.23  The  several  problems  of  propagation  discussed  in  this 
subsection  are  clearly-illustrated. 

For  each  release,  (Figure  3.4),  thermal  ionization  in  the  exploding 
thermite  charge  caused  immediate  signal  loss  followed  by  signal  recovery 
as  the  explosion  products  cooled  and  expanded.  The  first  release  was  at 
such  a low  altitude  that  the  barium  cloud  was  confined  to  a small  volume 
and  the  "footprint"  of  the  cloud  rather  quickly  moveu  away  from  the  ground 
receiver  due  to  rocket  motion.  Fhotoionization  of  the  barium  atoms  has 
a time  constant  of  approximately  15  seconds  and  the  record  shows  the  onset 
of  multipath  for  signals  reaching  the  receiver  a;  approximately  15  seconds 
after  the  release  and  continuing  for  approximately  30  seconds,  by  which 
time  the  receiver  viewed  the  direct  rays  only.  The  second  release  expanded 
to  a substantially  larger  volume  and  the  "footprint"  for  that  release  covered 
the  ground  receiver  for  approximately  35  seconds.  The  record  from  135-165 
seconds  illustrates  multipath  (higher  frequency  oscillations)  and  tht  signal 
attenuation  arising  from  defocusing  by  the  plasma  "ball".  The  increase 
in  attenuation  with  increasing  time  simply  reflects  the  photoionization 
time  constant  for  barium  atoms.  Shortly  after  165  seconds,  the  rocket  transmitter 
"comes  out  from  behind"  the  ion  cloud  and  the  receiver  observes  strong 
scintillation  arising  from  diffraction  and  multipath  effects.  At  the  same 
time,  the  phase-lock  tracking  loops  employed  in  the  receivers  lost  lock 
indicating  phase  changes  too  rapid  for  the  receivers  to  follow.  As  a rough 
approximation,  phase  rates-of-change  of  at  least  700  radians  per  second 
were  encountered. 

Experiments  similar  to  the  HAVEN  BARTER  experiment  just  described  were 
conducted  during  the  SECEDE  II  project  in  the  winter  of  1971  (see  Table  1.3). 

For  these  experiments,  a multi-frequency  transmitter  was  carried  by  a separate 
rocket,  and  occultat ions  by  the  barium  ion  clouds  occurred  substantially 
later  than  the  HAVEN  BARTER  occupations . Again  similar  effects  were  ob- 
served with  severe  phase  and  amplitude  scintillation  just  prior  to  and 
after  occupation.  For  the  SECEDE  II  releases,  four  receiving  stations  were 
operated  and  it  was  possible,  because  of  spatial  separations,  to  model  the 
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ion  cloud  as  seen  at  radio  frequencies.  The  model  was  constructed  by  assum- 
ing an  ellipsoidal  ion  cloud  distribution  with  self-consistent  cross  sections 
as  observed  from  the  four  receiving  ground  stations.24  The  geometrical  and 
ionization  parameters  of  the  ellipsoid  were  then  adjusted  to  match  the  dis- 
persive phase  measurements  for  each  of  the  stations.  The  resulLs  of  this 
modeling  procedure  are  shown  in  Figures  3.5  and  3.6.  The  EKL  experiments 
in  SECEDE  II  illustrate  a powerful  method  of  determining  the  cifects  of 
localized  plasma  concentration  and  distribution  on  radio  signal  transmissions. 
The  process  of  constructing  a radiofrequency  model  provides  a credible 


IV  EXAMPLE  SELECTED  FOR  ANALYSIS 


In  the  previous  section  the  several  kinds  of  propagation  problems 
giving  rise  to  communications  degradation  were  discussed.  Each  of  the 
problems  can  be  illustrated  in  the  propagation  experiments  conducted  by 
the  BRL;  the  nature  of  the  experiments  provided  data  to  permit  a quantitative 
estimate  of  effects  on  communications  systems  at  radiofrequencies  employed 
during  the  experiments  and  at  other  frequencies  as  well.  In  this  section, 
one  propagation  problem  is  analyzed  using  results  from  the  BRL  experiments. 

A second  problem  quite  appropriate  for  analysis  with  the  BRL  data  is  that  of 
collision  frequencies  in  the  lower  D-region.  Collision  frequencies  are  in 
themselves  not  a problem  but  their  proper  value  can  become  quite  important 
for  estimation  of  signal  absorption  in  the  presence  of  high  electron 
densities  in  the  lower  D-region.  Such  high  densities  did  occur  in  the  polar 
cap  absorption  (PCA)  event  in  November  196911  (and  see  Figure  4.4).  Using 
signal  absorption  measured  during  the  event,25  it  was  possible  to  compare 
calculated  and  measured  differential  absorption,  but  the  comparison  was 
carried  out  for  just  one  frequency  at  one  time  and  with  an  assumed  electron- 
neutral  collision  frequency.26  Other  work2 7 suggests  that  collision  fre- 
quencies could  vary  by  a factor  of  three  in  the  lower  D-region;  in  the  pre- 
sence of  high  electron  densities,  this  variation  can  markedly  alter  signal 
absorption.  A considerable  number  of  absorption  measurements  were  made 
by  the  BRL  during  PCA  69  and  these  data  probably  warrant  further  analysis 
in  the  context  of  looking  for  collision  frequency  variability  (if  such  exists) 
and  for  the  best  numerical  values  corresponding  to  conditions  during  PCA  69.* 

Focusing  and  Defocusing  of  Signals 

In  the  natural  ionosphere,  the  focusing  and  defocusing  of  radio  signals 
result  from  propagation  through  plasma  ''holes"  or  "blobs",  respectively. 

The  magnitude  of  the  effect  depends  primarily  on  ionization  gradients  normal 
to  the  direction  of  propagation.  In  the  natural  ionosphere,  the  strongest 

k 

Since  the  preparation  of  this  brief  report,  attention  has  been  called  to  a 
recent  report  by  W.  Swider  and  I.L.  Chidsey  in  which  the  radiofrequency 
absorption  data  were  analyzed  in  context  of  the  further  study  suggested  here. 
("Calculated  and  Measured  HF/VHF  Absorption  for  the  2-5  November  1969  Solar 
Proton  Event,"  Air  Force  Geophysics  Laboratory  TR-76-0053,  March  1976.) 
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gradients  most  likely  occur  in  the  auroral  region  where  strong  variations, 
in  time  and  place, occur  in  the  deposition  of  ionizing  particles. 

During  the  PCA  69  experiments  at  Ft.  Churchill,  Canada,  signal  focusing 
and  defocusing  was  encountered  on  two  occasions.  Data  from  the  propagation 
experiment  are  analyzed  in  this  section  to  provide  a semiquantitative  pic- 
ture of  postulated  plasma  blobs  in  the  auroral  ionosphere.  Such  a picture 
is  useful  only  to  indicate  the  ionization  gradients  possible  and  the  size  of 
the  ionization  inhomogeneities.  The  picture  does  provide  a model  for  esti- 
mating effects  on  communications  systems  operating  at  frequencies  other  than 
those  used  in  the  BRL  experiments  and  for  geometries  for  which  transmission 
through  a number  of  blobs  may  occur. 

As  indicated  earlier,  a ray-trace  analysis  is  required  for  accurate 
representation  of  signal  attenuation  under  conditions  of  defocusing  by  a 
complex  plasma  blob.  Such  an  analysis  is  beyond  the  scope  of  this  report 
in  the  context  of  deriving  a semiaccurate  picture  of  the  ionization  in- 
homogeneities encountered  in  PCA  69.  Rather,  an  approximate  analysis  method 
is  used  which  does  provide  a useful  insight  into  the  conditions  en- 
countered. This  simple  analysis  provides  answers  in  good  agreement  with 
measured  attenuations  for  the  HAVEN  BARTER  barium  release  and  for  which  a 
ray-trace  analysis  was  carried  out. 


In  the  HAVEN  BARTER  experiment,  the  combined  dispersive  doppler  and 
attenuation  measurements  gave  the  following  derived  parameters  for  the  barium 
plasma  cloud  (gaussian  distribution  of  Ba  ions  assumed):23 


and 


Time  constant,  x,  = 15  sec. 

23 

Limiting  Total  Cloud  Inventory,  TI ^ ^ = 7.0  x 10  ions  (or  electrons) 

2 2 

Scale  Radius  squared,  h = 4.3  + 0.031(t),  km" 

(j 


N - N exp\/h  /,  electron  distribution  in  cloud  as  tunction  of  radii 


h.2.=  (h2)  + 4Dt 
o 


T 1 = 


- TI  . (1-exp  yt)  = ion  inventory  as  function  ■'!  time 


1 


— ■ ■ ■ — - 
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N = peak  electron  density  at  center  of  cloud 
o 

t“  time  after  release  of  the  barium  cloud 

D=  ion  diffusion  coefficient  (approximately 
-3  2 -1 

7.5  x 10  km  sec  as  derived  from  the  experiment) 

h = radius  of  cloud  at  which  l=d  /e 
G o 

For  the  simplified  analysis,  the  barium  plasma  was  assumed  to  ix  a 
sphere  of  uniform  density  and  radius  h^,.  The  density  is  given  by, 

N"=  Ti  {(4/3)  (it  h^)}""1 

Using  these  simplifying  assumptions  and  Che  schematic  shown  in  Figure  4.1, 
an  expression  for  signal  attenuation  in  transmission  through  the  spherical 
plasma  can  be  derived.  Referring  tol’igure  4.1,  the  following  relationships 
hold: 

sin  a = L (sin  3)/R  (4.1) 

sin  [■  = sin  a e ^ (Snell's  Law)  (.4.2) 

- , 1/2 

p index  of  refraction  = [l-(80.5  N)/f“] 

-3 

N = electron  density  (m  ) 
f = transmission  frequency  (hz) 


s in  ,3  = 

sin  a 

[l-(80.5N/f2)]“1/2 

(4. 

2) 

- 

a - 6 

(4. 

3) 

Y = 

i 

o 

cc 

(4. 

4) 

d'  = 

L8U-  <- 

(Y+, ) = 0+2  (3-a) 

(4. 

5) 

AH  « 

R sin 

a/sin  d-* 

(4. 

b) 

Attenuation  ot  tm  radio  signal  transiting  the  spherical  plasma  is  calcu- 
lated by  noting  that  the  signal  power  radiated  i . . the  small  solid  angle, 
represented  by  , is  .pread  by  refraction  inti  a solid  angle  represented  by 
Reduction  in  signal  power  by  such  refraction  is,  simply, 
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P2/Pi  - A1^A2 


P2/P1  = | (L+H)2tan2e(  |(H+AH)2tan23‘*|_1 


(4.7) 


where  and  A„  are  areas  on  the  ground  illuminated  by  energy  contained  in 
the  solid  angles  represented  by  and  Q';  and  F.,  are  the  received 
power  per  unit  area  in  the  absence  and  presence  of  the  spherical  plasma, 
respectively.  The  signal  attenuation,  then,  is  given  by, 


Attenuation  (db)  = 10  log  j (L+H)  “tan^u  i j(H+AH)  “tan‘'0'*|  P (4.8) 


The  result  predicted  by  Equation  4.8  may  be  seriously  questioned  in 
light  of  the  radical  assumption  of  uniform  electron  density  distribution  in 
a refracting  sphere  rather  than  a quasi-gaussian  distribution  in  a plasma 
more  closely  resembling  a prolate  spheroid.  Some  compensation  does  occur 
in  that  the  radius  of  the  real  plasma  is  somewhat  greater  than  the  radius 
assumed  so  that  total  refraction  of  the  transiting  radio  signal  is  approxi- 
mately the  same  for  the  real  and  assumed  plasma  distribution.  Further,  the 

12  -3 

peak  electron  density,  N =9x10  m , derived  from  the  ray-trace  model  is 

° * 13  -3 

reasonably  close  to  the  N =1.3x10  m used  in  the  approximation. 


Tlie  attenuation  predicted  by  Equation  4.8  is  compared  to  measured 
data  and  the  ray-trace  analyses  in  Table  4.1.  Measured  data,  comparable 
to  those  shown  in  Figure  3.4,  ire  available  for  the  36  and  146  MHz  trans- 
missions.23 The  comparison  for  146  MHz  is  not  valid  in  that  measured 
attenuation  was  effectively  zero  at  times  corresponding  to  maximum  attenua- 
tion at  36  and  73  MHz,  The  maximum  measured  attenuation  at  146  MHz  occurred 
at  approximately  10  seconds  after  plasma  release  at  which  time  measured 
attenuations  at  36  and  73  MHz  were  approximately  half  their  maximum  value. 

In  effect,  the  "footprint”  of  the  plasma  cloud  for  the  146  MHz  transmission 
did  not  cover  the  receiving  station  at  30  seconds  after  plasma  release. 

An  examination  of  attenuation  predicted  by  Equation  4.8,  and  with  refer- 
ence to  Figure  4.1,  shows  that  the  most  sensitive  parameter,  by  far,  is  the 
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Table  4.1 


MEASURED  vs  PREDICTED  MAXIMUM 
(Approximately 

ATTENUATION,  HAVEN  BARTER  RELEASE  No, 
30  Seconds  After  Release) 

. 2 

Frequency1 

36  MHz 

73  MHz 

146  MHz 

Measured 

-30+?r,  ,!b 
-10 

-12  + 2 db 

-2.5  ± 1 db" 

Ray-Trace  Analysis 

-28  ± 1.2  db 

-8.3  f O.u  db 

— 

Predicted  by  Equation 
4.8 

-33  db 

-12  db 

-5  db 

1 Rounded  to  nearest  MHz. 


Measured  at  approximately  10  seconds  after  release. 
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plasma  electron  density.  For  the  realistic  case  in  which  a uniform  electron 
density  is  replaced  by  a plasma  with  density  gradients,  the  ionization 
gradient  becomes  an  important  parameter.  Removing  the  transmitter  to  greater 
distances  from  the  plasma  cloud  (parameter  L in  Figure  4.1)  is  effective 
in  reducing  the  "footprint"  of  the  cloud  as  seen  from  the  ground.  The 
effect  of  increased  L in  Equation  4.8  is  compensated  by  a decrease  in 
so  that  attenuation  within  the  footprint  remains  approximately  the  same. 

Attenuation  data  representing  plasma  defocusing  for  HAVEN'  BARTER, 

SECEDE  III  and  PCA  69  are  shown  in  Figure  4.2.‘3  11  10  The  attenuation  valu> 
shown  for  HAVEN  BARTER  and  the  GUM  and  ELM-2  releases  in  SECEDE  111  ae 
reasonably  consistent  in  that  larger  barium  payloads  were  used  in  Lh<. 

SECEDE  experiments.  The  attenuation  for  ELM-1  appears  inconsistent  in  that 
the  barium  ion  yields  for  ELM-1  and  ELM-2  were  very  similar.  The  most 
likely  explanation  for  the  marked  differences  in  the  ELM  releases  derives 
from  the  different  release  altitudes.  EIM-1  was  released  at  140  km  anc 
confinement  of  the  release  by  the  ambient  atmosphere  was  substantially 
greater  than  for  the  higher  release.  Even  so,  the  maximum  electron  densit\ 
indicated  for  ELM-1  was  close  tc  10  m ' , a value  perhaps  thret  time  greater 
than  peak  densities  (mostly  inferred)  in  other  barium  releases. 

For  a natural  disturbance  in  the  auroral  regions,  sucli  as  that  produced 
^ by  a polar  cap  absorption  (PCA)  event,  the  simple  picture  of  a "plasma  ball" 

giving  rise  to  abnormal  attenuation  is  scarcely  adequate.  Figure  4.3 11  I 

shows  measured  attenuation  during  flight  of  an  early  rocket  (B-2)  in  th< 

PCA  event  of  November  1969.  The  absence  of  appreciable  attenuation  befot 
125  seconds  of  flight  indicates  a relatively  low  concentration  of  electron- 
below  80  kilometers.  Electron  densities  derived  from  dispersive  phase 
measurements  are  shown  in  Figure  4.4  and  substantiate  the  low  attenuatin' 

(absorption  in  this  case)  measured  on  flight  B-2.  Beginning  at  ap, roximatel > 

125  seconds  (tj)  a marked  enhancement  (focusing)  of  the  signal  occur  ^ 
followed  at  140  seconds  (t, ) by  the  onset  ot  strong  attenuation  (defe'eusing) 
of  the  signal.  The  period  of  abnormal  attenuation  may  be  interpreted  as 
rocket  overflight  of  a small  region  of  enhanced  electron  density  arising 
1 rom  particli  deposition  below  the  rocket.  A highly  schematic  repiesentat-  n 
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Figure  4.2  ATTENUATION  FROM  SIGNAL  DEFOCUSING 
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Figure  4.3  RADIO  FREQUENCY  ATTENUATION  ILLUSTRATING  FOCUS/DEFOCUS 


of  the  disturbance  cross  section  (in  the  trajectory  plane)  is  shown  in 
Figures  4.5  and  4.6.  The  times  t^,  t^,  etc.,  shown  in  Figures  4.3  and  4.6 
roughly  correspond.  The  height  of  the  disturbance,  h^,  is  constrained  in 
that  it  must  be  below  the  rocket  but  above  the  region  (above  80  km)  of 
large  absorption;  the  size  is  constrained  by  the  short  duration  of  the  large 
signal  attenuation.  If,  for  example,  the  discontinuity  ("blob")  were  at  an 
altitude  of  90  km,  the  horizontal  motion  of  a raypath  at  that  altitude  is 
appro ximately  2.5  km  between  120  and  160  seconds  of  flight  time.  There  is 
a major  constraint,  also,  on  the  electron  density  in  the  discontinuity  in 
view  of  the  data  shown  in  Figure  4.4.  These  several  constraints  make  it 
possible  to  estimate  the  size  and  character  of  the  discontinuity  despite 
the  highly  fictitious  representation  shown  in  Figure  4.1  and  the  attenua- 
tion expression  derived  therefrom. 

2 2 , 

From  Equation  4.8,  the  ratio  (tan  0)/(tan  0 ) clearly  dominates  in 
determining  the  attenuation.  From  Equations  4.1,  4.2,  4.5,  4.8 
and  the  measured  attenuation  of  -26  db  (Figure  4.3)  for  36  MHz,  a "blob" 
density  of  3-5x10^  electrons  m ^ and  a "blob"  diameter  of  200  meters  at 
10-15  kilometers  below  the  rocket  provides  an  approximate  fit  to  tin  data. 

The  calculated  attenuation  for  this  representation  of  the  discontinuity  is 
extremely  sensitive  to  electron  densities  assumed  and  the  size  of  the  dis- 
continuity. Larger  electron  densities  are  essentially  prohibited  by  the 
very  high  ionization  rates  required;  a density  of  3x10^  requires  highly 
localized  ionization  rates  of  the  order  of  10^*  ion  pairs  m * sec  1 , a 
value  approximately  5 times  greater  than  measured  by  other  experiments  during 
the  PCA  69  event.12  At  altitudes  below  90  km,  the  increasing  recombination 
rates  would  require  still  higher  ionization  rates  for  a given  electron  density. * 9 

The  thrust  oi  the  above  analytical  argument  and  representation  of  the 
ionization  discontinuity  is  not  that  the  measured  attenuations  an  provide  a 
picLur'.  of  the  ionization  structure  during  the  highly  disturbed  conditions 
encountered  in  the  auroral  regions  but  rather  that  there  is  a physical 
basis  for  the  observations.  From  the  standpoint  of  communication  the 
existence  of  blobby  structures  provides  a mechanism  for  signai  loss  in 
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addition  to  the  more  commonly  observed  absorption  of  signal  energy  arising 
from  relatively  high  ionization  at  lower  altitudes  and  attendent  electron- 
neutral  collisions.  Considering  the  fact  that  anomalous  attenuation  was 
observed  on  two  rocket  flights  at  Ft.  Churchill  and  probably  occurred  on 
one  of  the  flights  during  the  1969  SECEDE  III  experiments  at  Poker  Flat, 
Alaska,10  the  occurrence  of  signal  defocuslng  in  auroral  regions  would 
appear  statistically  significant. 

The  presence  of  sharp  discontinuities  In  Ionization  gives  rise  to 
another  potential  problem  for  phase  coherent  communications  systems.  For 
the  HAVEN  BARTER  plasma  releases,  the  phase-lock  loop  operating  at  36  MHz 
could  not  track  at  the  time  the  radio  signals  passed  the  edge  of  the  barium 
ion  cloud  (approximately  170  seconds  in  Figure  3.4).  This  loss  of  phase 
coherence  arises  mainly  from  multipath  for  signals  grazing  the  region  of 
high  ionization.  A similar  situation  arose  during  flight  of  the  B-2  rocket 
in  the  PCA  69  experiments.  For  this  flight,  dispersive  phase  was  obtained 
from  the  coherent  pairs  of  frequencies,  36/145  MHz  and  145/583  MHz.  Starting 
at  approximately  130  seconds  (see  Figure  4.3),  phase  lock  was  lost  for  the 
36  MHz  signal,  most  likely  from  seve.o  multipath,  with  an  accumulated  phase 
loss  of  approximately  30  radians  which  was  never  recovered. 

Certainly  the  most  spectacular  example  of  a highly  disturbed  natural 
ionosphere  during  any  period  of  observations  by  the  BRL  group  occurred 
during  the  IGY.  At  that  time,  ionospheric  effects  were  determined  by  the 
relatively  crude  trajectory  comparison  method  and  Faraday  rotation  of  the 
two-way  DOVAP  signals.  The  DOVAP  system  was  tracking  an  Aerobee-Hi 
rocket  launched  at  0025  CST  from  Ft.  Churchill  on  24  November  1958.  The 
vehicle  flew  into  a weak  aurora  but  the  effect  on  thi  tracking  system  was 
dramatic.  Rapid  phase  shifts  and  fading  of  the  37  MHz  signal.-,  degraded 
tracking  accuracies  to  the  extent  that  the  trajectory  comparison  method 
could  not  be  used.  Faraday  rotation  data  indicated  onset  and  overflight  ot 

regions  of  extremely  enhanced  electron  density  (both  up-leg  and  down-leg) 

. , 13-3 

with  very  sharp  density  gradients.  An  estimate  of  2x10  electrons  v 

in  a very  small  volume,  and  at  an  altitude  less  than  120  kni,  give  son,  mat  h 

to  the  observed  data.18 
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Knowledge  of  the  existence  and  some  quantitative  measure  of  attenua- 
tion and  multipath  associated  with  blobby  plasma  structures  in  the  natural 
ionosphere  makes  possible  compensations  in  communications  systems  operating 
in  that  environment.  Much  more  difficult  and  potentially  serious  are  these 
same  effects  on  signal  strength  and  signal  phase  in  the  severe  environment 
produced  by  high  altitude  nuclear  events. 

I 

[- 
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V SUMMARY  AND  CONCLUSIONS 


Previous  sections  have  given  a most  cursory  review  of  efforts  at  the 
Ballistic  Research  Laboratories  in  their  studies  of  radio  propagation 
through  the  ionosphere.  Following  rocket  experiments  in  Kauai  in  1971, 
personnel  and  work  reassignments  together  with  a disappearing  mission  led 
to  termination  of  upper  atmosphere  research  with  rocket  sounding  techniques. 
A vigorous  laboratory  program  for  studying  atmospheric  chemistry  reactions 
relevant  to  altitudes  of  interest  to  the  BHD  program  was  continued  and 
complex  balloon-borne  experiments  were  initiated.  However,  a capability 
for  atmospheric  research  with  rockets  by  the  BRL  is  probably  beyond  recall. 

There  have  been  many  contributions  by  the  BRL  in  propagation  research 
not  touched  upon  or  referenced  in  this  report.  Work  in  antenna  theory  and 
design  together  with  the  development  of  equipments  for  the  reception  and 
processing  of  extremely  weak  signals  was  found  in  the  vanguard  of  these 
technologies.  The  development  of  digital  ray  tracing  techniques  was  an 
important  contribution.  The  use  of  dispersive  doppler  information  from 
spaced  ground  receivers  to  describe  plasma  structure  was  an  unusual  applica- 
tion of  the  BRL  propagation  studies. 

In  the  long  run,  BRL  measurements  of  propagation  effects  in  the  auroral 
and  nuclear  disturbed  ionospheres  will  constitute  an  important  if  not 
major  part  of  the  available  data  base  for  such  disturbances.  Nuclear  tests 
in  the  atmosphere  are  not  foreseeable  in  the  near  future,  and  there  is  no 
assurance  that  natural  disturbances  of  the  magnitude  seen  in  IGY  will  be 
available  against  which  to  test  communications  systems  of  the  future. 
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